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ABSTRACT 

Context. Imaging polarimetry is a powerful tool for detecting and characterizing exoplanets and circumstellar environments. 
Polarimetry allows a separation of the light coming from an unpolarized source such as a star and the polarized source such as a 
planet or a protoplanetary disk. Future facilities like SPHERE at the VLT or EPICS at the E-ELT will incorporate imaging polarime- 
try to detect exoplanets. The Extreme Polarimeter (ExPo) is a dual-beam imaging polarimeter that currently can reach contrast ratios 
of 10', enough to characterize circumstellar environments. 

Aims. We present the data reduction steps for a dual-beam imaging polarimeter that can reach contrast ratios of 10^. 

Methods. The data obtained with ExPo at the William Herschel Telescope (WHT) are analyzed. Instrumental artifacts and noise 

sources are discussed for an unpolarized star and for a protoplanetary disk (AB Aurigae). 

Results. The combination of fast modulation and dual-beam techniques allow us to minimize instrumental artifacts. A proper data 
processing and alignment of the images is fundamental when dealing with large contrasts. Imaging polarimetry proves to be a powerful 
method to resolve circumstellar environments even without a coronagraph mask or an Adaptive Optics system. 

Key words. Instrumentation: polarimeters - Polarization - circumstellar matter 



1. Introduction 

The direct detection and characterization of exoplanets is one of 
the main goals of the next generation of instruments and tele- 
scopes, such as the E-ELT ( |Hook|2009) l. Despite the increasing 
amount of detected exoplanets (e.g., http://www.exoplanet.eu/), 
only a few of them have been directly imaged ([Marois et al. 
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2008) Kalas et al. 2008 Lagrange et al. 2010 Todorov et al. 



Polarization is a powerful tool for detecting and character- 
izing circumstellar environments, such as protoplanetary disks, 
debris disks and exoplanets (Kuhn et al. 2001 ; Stam et al. 2005 , 
|2006| l. While the light coming from the star is largely unpolar- 
ized, the scattered light from exoplanets and circumstellar envi- 
ronments is highly (linearly) polarized. This makes it easy to 
separate the planet signal from the unpolarized starlight, and 
reach the required, large contrast ratios. 



et al 



Future exoplanet instruments such as SPHERE (Thalmann 
2008| l, GPI ( [Graham et all2007l ) or EPICS ( Kasper et al 



2010 will include imaging polarimetry, spectropolarimetry, or 
both, to characterize exoplanets. To reach the high contrast ra- 
tios that are necessary to detect an exoplanet, these instruments 
will combine coronagraphs. Adaptive Optics and polarimetry. 
The combination of the first s two elements ca n reduce the stellar 
light by a factor of 10^ (Martinez et al. 2010|, while polarimetry 



alone can reach contrast ratios of 10^ (Rodenhuis et al. ] |2011 



The identification of instrumental artifacts and spurious signals 
that appear at these extremely high ratios is fundamental to un- 
derstand the data produced by the new generation of exoplanet 
imagers. 

ExPo is a dual-beam imaging polarimeter that combines 
speckle imaging with dual-beam techniques, currently reaching 



contrast ratios of 10^ without a coronagraph or adaptive optics 
system (Rodenhuis et al. 2008 201 1 1. The sensitivity reached by 
this instrument allows us to image circumstellar environments 
such as protoplanetary disks or dust shells with unprecedented 
accuracy (e.g. [Canovas et al.|201 la|b| l 

In this paper we present the data reduction approach for high 
dynamic range polarimetry. The whole data-reduction process 
is shown in Fig. [T] Results show that the dual-beam technique 
effectively minimizes systematic errors, leading to sensitive po- 
larization measurements of different circumstellar environments. 

In the following sections, the instrument is briefly described, 
followed by the analysis of the instrumental artifacts. The image 
reduction process and a discussion of first results are presented 
at the end. Preliminary results of the well-known protoplanetary 
disk AB Aurigae and the diskless star HD12815 are shown. 



2. Instrument Description 

ExPo is a dual-beam imaging polarimeter that currently reaches 
contrast ratios of 10^. It is designed to characterize circumstel- 
lar environments, such as the disks around Herbig Ae/Be and T 
Tauri stars. ExPo is currently a visitor instrument at the William 
Herschel Telescope (WHT), where it has been used for four suc- 
cessful observing campaigns. 

ExPo has a wavelength range from 400-900 nm, and its field 
of view is 20" x 20". It works with exposure times of 0.028 sec- 
onds to minimize the effects of atmospheric seeing. Systematic 
errors are minimized by modulating two opposite polarization 
states. The combination of dual-beam, short exposures times and 
polarization modulation allows us to reach contrast ratios of 10^ 
without the aid of Adaptive Optics (AO) or coronagraphs. 
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As a dual-beam instrument, the incoming beam of light is 
divided into two orthogonally linearly polarized beams. As a re- 
sult of this, two simultaneous images with opposite polarization 
states can be recorded on the CCD at the same time. 

The schematic layout of the instrument is shown in Fig. |2] 
The light beam enters the instrument from the left. The disper- 
sion introduced by the atmosphere (e.g.'Wynne'1993 ) is reduced 
by an Atmospheric Dispersion Corrector (ADC). The field stop 
prevents beam overlap on the detector. The filter wheel contains 
broadband and narrowband filters as well as an opaque filter to 
obtain dark measurements. The polarization compensator mini- 
mizes the instrumental polarization introduced by the telescope. 
The light beam is currently modulated by a single Ferroelectric 
Liquid Crystal (FLC), which switches between two states (A and 
B) with orthogonal polarization states. Future versions of this in- 
strument may use an achromatic 3-FLCs Pancharatnam modula- 
tor ( |Pancharatnam 1 955 ; Keller 2002 ), but this will not influence 
the data reduction approach described in here. The beamsplitter 
divides the incoming light into two beams (left and right) with 
perpendicular polarization states. These two beams are finally 
imaged onto the EM-CCD camera, which records two simul- 
taneous images with orthogonal polarization states. During the 
observation, the whole instrument remains static, with only the 
FLC switching from one state to the other. 



2.1. ExPo Images 

With the FLC in two different states, two different frames (A 
and B) are recorded. A total of four images (i.e. two images per 
frame, corresponding to the left and right beams) are produced 
for each FLC cycle: A^, A«, Bi, Bjf, where the subscripts L and R 
stand for left and right, respectively. Due to the FLC modulation, 
A/, and B« contain the same polarization information, as happens 
to A« and Bi. To account for all the components of these images, 
we use the following notation: 



FLC cycle can be described by: 



Ta,b 
Tl,r 



FLC transmission 
Beamsplitter transmission 



f tofAL,AR,BL,BR Fl^t Fields 



PSF 

N 



Intensity, Polarization images 
Point Spread Function 
Noise 



where Ta,b accounts for the FLC transmission in its two pos- 
sible states, and Tl,r accounts for any transmission difference 
between the two beams. /^/flfAL,AR,BL,BR accounts for the differ- 
ent flat fields on the left and right sides of the CCD for the A and 
B images. The noise term includes the read-out and the photon 
noise. The polarization image includes any polarized light that 
enters the instrument, including sky polarization and polariza- 
tion introduced by the telescope. The PSF can be expressed as 
the convolution of the "optics" and (PSFopt) and "seeing" (S) 
terms of die PSF (e.g. |Roddier|2004[ ): 



PSF = PSFopt * S, 



(1) 



where PSFopt accounts for the static optical elements of the in- 
strument and telescope, and S accounts for the dynamic atmo- 
sphere. Any small imperfection in the beamsplitter, or difference 
in the optical path between the two beams will translate into dif- 
ferent PSFopt for the left and right images: PSFl and PSFr. On 
the other hand, the atmosphere will change with time, so S will 
also change for the A and B frames: S a and Sb- Taking all of 
this into account, the four images recorded by the CCD in one 





= Ta 


Tl 


FlatAL 


0.5 


Ar 


= Ta 


Tr 


FlatAR 


■0.5 


Bl 


= Tb 


Tl 


FlatBL 


0.5 


Br 


= Tb 


Tr 


FlatBR 


0.5 



where Al stands for the left image taken when the FLC is in its 
A state, and so on. Nal refers to the noise associated with this 
image. Symbols "■" and "*" stand for multiplication and convo- 
lution, respectively. The factor 0.5 comes from the beamsplitter: 
the total incoming beam is divided into two. A "pure" intensity 
image can then be generated by adding the left and right beams. 

Once a data set is recorded, the FLC is rotated by 22.5°, and 
a new data set is acquired. Since ExPo does not have a derota- 
tor, the field rotation will not only rotate the image, but also the 
polarization plane of the incoming light. To avoid this, the im- 
ages are calibrated according to the reference system of the ob- 
served target, and then transformed to the reference system of the 
instrument. The Stokes parameters Q and U are obtained after 
combining and calibrating ( Rodenhuis et al.|20l 1 ) four datasets 
with flie FLC rotated by 0°, 22.5°, 45° and 67.5° (e.g. |Patat &| 



Romaniellol (2006,) ). Two FLC positions (0° and 22.5° or 45° 
and 67.5°) are necessary to produce calibrated images, so two 
redundant data sets are obtained after calibrating. Random er- 
rors (e.g. uncorrected cosmic rays, FLC ghosts) are removed by 
comparing these two data sets. 

3. Image Preparation 

In this section we describe the initial reduction steps that must 
be applied to the raw data. Once the images are corrected for 
dark, flat and smearing effects, they are ready to be aligned and 
combined as explained in §4 and §5. 

3.1. Dark Reduction 

A dark image is an image taken with a given exposure time and 
with the CCD shutter closed. However, this is not possible due to 
the short exposure times required by ExPo, as will be explained 
in the next subsection, and the CCD shutter must remain open 
during the whole observation. An opaque filter is placed in the 
filter wheel to measure dark frames. 

Every dark frame is checked to detect and remove Cosmic 
Rays (CR). CR are identified as pixels showing values greater 
than a threshold value. This threshold is defined according to the 
characteristics of the dark observation, i.e. dark frames obtained 
with different CCD gain values will require different thresh- 
old values to detect the CR events. This process produces good 
results when reducing dark frames. Another approach is used 
when reducing the final science images, as is explained in §3.4. 
Finally, a very low-noise master dark is generated by combining 
at least 8000 dark frames (equivalent to a total integration time 
of roughly 4 minutes). 

3.2. Smearing and Second-order Corrections 

Due to the short exposure times used by ExPo, the CCD runs 
in frame-transfer mode, which means that the detector is read- 
ing and recording frames at the same time. To achieve this, the 
CCD chip is split in half: one half for exposing and another half 
for reading. Once an image is recorded on the "exposing" half. 
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Fig. 1. Pipeline Block Diagram. 
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Fig. 2. Schematic layout of the instrument in the Nasmyth focus of a telescope. 



is shifted to the read-out area. The shutter remains open during 
the whole process, and this, unfortunately, introduces spurious 
effects that must be removed in the initial reduction of the raw 
images. 

Each time a frame is transferred to the read-out area of the 
CCD, there is some leftover charge. As a result of this, an extra 
background is added to the next recorded frame. Moreover, in 
case the observed target is too bright, a streak in the direction of 
the charge shifting will appear when exposing for a long time, or 
averaging many short-exposure images. 

Due to the field mask, some areas of the CCD are not illu- 
minated (Fig.|3]l. Any signal measured in these non-illuminated 
regions must come either from the extra background mentioned 
above, or from the read-out noise. The information contained 
in these "dark" pixels is used to remove these two effects. An 



example of this is shown in Fig. |4] the spurious signal is first 
measured outside of the field mask, and then subtracted from the 
raw image. 



3.3. Flat Reduction 

A set of dome flats is taken at the beginning and the end of each 
observing night. A minimum number of 100 frames are recorded 
for each of the four FLC positions (0°, 22.5°, 45°, 67.5°). Eveiy 
frame is then corrected for dark and smearing effects. As with the 
images, it is convenient to distinguish between Flat^^, Flatp^^, 
FlatBL and FlatsR. Any polarization introduced by the dome 
lights, beamsplitter and FLC transmission differences disappears 
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Fig. 3. Normalized Flat Fields with two different field masks. 
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Fig. 4. Left: Raw image, plotted at very high contrast. Two spu- 
rious ghosts can be seen above the star (discussed later in §6.4) 
The horizontal white line represents the plot that is shown on the 
right. Right: The black line shows the horizontal cut before the 
correction, the green line shows the smearing correction effect 
described in §3.2. 



when normalizing each of these flats independently. All of the 
flat field images are combined to generate four master flats. 

The analysis of the flat fields shows that the CCD has no 
dead pixels. 

3.4. Cosmic Ray Correction 

Before the images are combined and analyzed, they are corrected 
for cosmic rays (CR). A standard ExPo data set comprises 20000 
images (around 10 minutes exposure time) per FLC position. 
The CR contribution cannot be neglected since we are dealing 
with contrast ratios of a; 10^. Due to the nature of the ExPo im- 
ages, with very short exposure times and a FOV of 20" x 20", 
it is not trivial to discriminate between a real speckle and a CR. 
Standard CR rejection methods are based on comparing an av- 
erage image with the individual ones, and then apply a sigma- 
rejection algorithm. However, the speckle pattern disappears in 
the average; the average and individual images are indeed quite 
different. A different approach, based on the advantage of having 
two simultaneous images in the same frame, is used in this anal- 
ysis. A box of 5 X 5 pixels is extracted around the brightest pixel 
of the left and right images. The standard deviations (cr) of these 
two boxes are then compared. Experimental results show that, 
if the brightest pixel corresponds to a CR and not to a speckle, 
its associated cr will be at least 3 times greater than the cr of a 
true speckle. CR's are detected and corrected according to this 
threshold: if one of the boxes contains a CR, all its pixels are set 
to the mean value of the whole image. 

In case both the left and right images are affected by CR's, 
this method will still work because the probability that the cos- 
mic rays have the same standard deviation is very low. Fainter 
CR can be divided into two categories: the ones that fall inside 
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Fig. 5. Polarized images of AB Aurigae (seeing was ~ 1.4"). 
Setl shows a calibrated image generated with the FLC placed at 
0° and 22.5°. Set 2 shows the polarized image generated after 
placing the FLC at 45° and 67.5°. The white circle shows the 
position of a cosmic ray that does not appear in the first image. 



the speckle pattern and the ones that fall outside of it. The former 
are, in practice, impossible to discriminate from true speckles. 
The latter are removed from the final images at the end of the 
analysis. Figure [3] shows an example of this process: one uncor- 
rected CR appears in one set but not on the other, so it can be 
easily identified and removed from the final image. 



4. Image Alignment 

Once these corrections are applied, the images are aligned before 
they are combined. During the image acquisition, the telescope 
guides on the observed target (i.e. there is not any reference star 
to be tracked). As a result of this, the accuracy of the guiding 
decreases, specially in the fainter (my > 9) targets. Because of 
the small field of view of ExPo, any small telescope misalign- 
ment can shift the image center by several pixels. Additionally, 
due to the short exposures used, the PSF is a speckle pattern. All 
of these contributions make image alignment one of the main 
sources of error in the processing of ExPo data. To reduce mis- 
alignment effects, two centering algorithms were tested: bright- 
est speckle and template cross-correlation. 

The brightest speckle method (or shift-and-add method) is 
used in lucky imaging (e.g. Law et aL]|2006| l. The images are 
aligned according to the position of the brightest speckle in every 
image. 

The template cross-correlation method first requires a start- 
ing template or reference image. Once the template is provided, 
each single image is cross-correlated with it. The image center 
is now defined as the position of the maximum value of this op- 
eration. Different templates such as a simulated gaussian-shaped 
Point Spread Function (PSF) or an average speckle-centered real 
PSF were used in these calculations. 

Once the images are centered, and before they are combined, 
they are again aligned to minimize any other source of misalign- 
ment, such as sub-pixel or flat-field gradients. These effects are 
corrected with a custom version of the drizzle code of lFruchterl 
|& Hook| ( |2002 i In this code, each pixel is expanded into several 
pixels, making it possible to shift and align images with sub- 
pixel accuracy. 

The analysis shows that the best results are obtained when 
combining these three methods. First, a template is generated 
by centering each left and right image according to the brightest 
speckle. An average non-gaussian shaped PSF is generated after 
combining all these images. Then, this PSF is used as the input 
template for the second method: each individual image is now 
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Fig. 6. Polarized (uncalibrated) results for a diskless (unpolar- 
ized) star: HD 1228 15. Four different alignment approaches are 
shown. Left column: images are aligned with a template. Right 
column: images re aligned according to the brightest speckle 
(shift-and-add). Top row: No sub-pixel alignment. Bottom row: 
sub-pixel alignment. Color bar: amount of (uncalibrated) polar- 
ization, P. 

cross-correlated with this template. Finally, each 'right' image is 
aligned with respect to the 'left' one, by using the drizzle code. 
Each pixel is subdivided into 9 sub-pixels, which results in an 
alignment accuracy of a third of a pixel. Any shift introduced by 
the FLC is also corrected by this centering process. 

To check how misalignment affects the data reduction, sev- 
eral unpolarized (diskless) stars were analyzed. In these cases, 
no polarization pattern is expected to be measured. The accu- 
racy of the alignment can be quantified in terms of the stan- 
dard deviation cr of this pattern. Figure |6] shows the polariza- 
tion pattern measured for the unpolarized star HD122815. The 
top row shows the difference between the template-center and 
the speckle-center. The bottom row shows the same result, but 
now including sub-pixel alignment with a third of a pixel accu- 
racy. The combination of template center plus sub-pixel align- 
ment shows almost no structure at all (Fig.|6] lower-left corner). 

From Fig. |6] it is clear that most of the misalignment 
disappears when aligning the images according to the cross- 
correlation method. No significant improvement is detectable 
when aligning with sub-pixel accuracy due to the lack of struc- 
ture in a diskless star. 



5. Image Combination 



As described by Bagnulo et al.| ( 2009| l there are two methods 
of combining images produced with a dual-beam polarimeter to 
produce polarization images: one based on the image ratios and 
the other based on the image differences. For ExPo, these two 
approaches can be expressed in terms of Eq. |2]as 

P = 0.25 ■ f^?^ - l) ■ 0.5 ■ (Al + Ar + B^ + Br) (3) 

for the double-ratio method, and as 

P = 0.5 (Pa - Pb) = 0.5 ((Al - Ar) - (Bl - Br)) (4) 




-10 • -10 

[Arc SccsJ 

Fig. 7. AB Aurigae images, obtained after applying the Double 
Difference (left) and Double Ratio (right). Color Bar: amount of 
(uncalibrated) polarization, P. 



Notation 


Expression 


AA 


Al - Ar 


AB 


Bl - Sr 


AL 


Al-^l 


AR 


Ar-Br 


6i 


Tl ■ PSFl - Tr • PSFr 


<5p 


Tl ■ PSFl + Tr • PSFr 


s[ 


Ta ■ 5' A - Tb ■ 5 B 




Ta ■ 5 A + Tb ■ 5 B 


ANa 


A'al - A^AR 


ATVb 


A'bl - A^BR 



Table 1. Notation used during the development of the image 
analysis. 



for the double-difference method. Both methods were investi- 
gated. The results prove that better results are obtained when ap- 
plying the double-difference method instead of the double-ratio 
explained in detail in Keller dT996jl and Donati et al. (1997i. 
This result is shown in Fig. |7j where the same data set was 
reduced applying the double-differences method (left) and the 
double-ratio method (right). This can be explained by the dif- 
ference of the speckle pattern for each image. Even though Al 
and Ar are recorded simultaneously, they show small differences 
in their speckle patterns due to instrumental effects. Due to the 
huge dynamic range present in our images, the slightest differ- 
ence between two images will produce a big effect when calcu- 
lating a ratio. However, the difference method is less sensitive to 
this. Therefore, the latter method was chosen for ExPo. At this 
point the dual-beam + beam-exchange technique shows its ad- 
vantages. The combination of the four images generated in each 
FLC cycle minimizes the systematic errors. Following the nota- 
tion introduced in §2, all terms that appear when developing the 
equations are listed in Table[T] Once the images are corrected for 
dark, flat and smearing effects, the following differences can be 
defined: 



AA = 0.5 ■ Ta ■ {/ * 5i + P * 5p) * 5a + AA^A, 
AB = 0.5 ■ Tb ■ {/ * 5i - P * 6y] * 5b + AA^b, 



(5) 
(6) 



where 6i and 5p refer to differences between the left and right 
beams. The difference between Tl and Tr produces a bias be- 
tween left and right images, while the difference between PSFl 
and PSFr differentially distorts the images. In case the beam- 
splitter does not produce two identical beams (i.e. due to optical 
imperfections) 6\ will be non-zero. Since ExPo looks for polar- 
ized intensities of P a: 10"^'/, this term should be on the order of 
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Unpolarized Star 



AB Aurigae 




[Arc ■ SecsJ 

Fig. 8. Normalized difference between left and right beams. The 
asymmetrical pattern around the center means shape differences 
(i.e. PSFl and PSFr produce different patterns). The color bar is 
in units of 10"^, meaning that this is the magnitude of the beam 
differences. 



« 10"^^ or smaller, as will be explained in the next paragraphs. 
As soon as this condition is not satisfied, the polarization values 
will be affected by the uncorrected intensity differences. 

To evaluate the effect of 6i on the polarimetry results, two 
average intensity "left" and "right" images were produced by 
averaging left and right beams: /l = Al + Bl.^r = + ^r- 
Both images were aligned with a third of a pixel accuracy, then 
normalized and finally subtracted. The noise contribution can 
be neglected in this calculation, since its contribution is several 
orders of magnitude below the intensity level. The result is de- 
scribed by: 



II - h - (I * + P * 6^) * 6i, 



(7) 



where 6^ and dp refer to any difference between the A and B 
frames. The difference between Ta and Tb will affect the whole 
image like a bias value, and the difference between S a and 
will affect the image shape. The later will be averaged out when 
combining enough images, i.e. the seeing effect will affect the 
A and B frames in the same way. FLC transmi ssion differences 
are n ot higher than 10"^ for our current FLC (Rodenhuis et al. 
201 1| , therefore 6[ ^ 10 ^ and will be very close to unity 



dp 1 . Taking all of this into account, the first term of Eq. 7 
is much larger than the second term: I * 6p » P * (J^, so Eq. 
becomes: 



(8) 



Figure |8] illustrates this result for the diskless star HD 1228 15, 
observed with the FLC oriented at 0°. The image is scaled to 
its maximum, i.e. the beam differences have an amplitude of 
=s +4- 10""^, which means that this is the magnitude of Sj. Similar 
values are found when repeating this experiment for different 
FLC orientations. As can be clearly seen from these results, 
beamsplitter imperfections are an important effect that needs to 
be taken into account when dealing with contrast ratios of 10"* or 
higher. Any differences in the PSF of each beam can modify the 




[Arc SccsJ 

Fig. 10. Polarization Degree for the unpolarized (diskless) star 
HD122815 (left), and for AB Aurigae (right). 



final result when aiming for such high contrast ratios. This con- 
tribution is, however, minimized when working with the double- 
difference method: 

n njl 

^(AA-AB) = ^ (/ * 5i * 5; -H P * 5p * (5p -I- AA^A - AA^b) -(9) 

In this case, the four images produced in each FLC cycle are 
first combined, and then averaged over the whole data set. Now 
the intensity term is convolved with the convolution of 61 and 
^j. The latter will have a magnitude of 61 * 6^ ~ 10"^, so the 
intensity will be decreased by this factor Both 5p and 5'^ are 
close to unity, so their convolution will be also close to unity. The 
polarization term is then not increased nor decreased. Finally, the 
noise term is not modulated by any delta factor Its contribution 
is minimized by the double-difference, and it will be averaged 
out. This term will define the maximum contrast attainable with 
ExPo. 

Figure |9] summarizes the previous results for both a diskless 
star (HD122815) and a protoplanetary disk (AB Aurigae). The 
first two columns starting from the left show the result given by 
Eq.|5]and Eq.|6] respectively, when applied to these two targets. 
The result of a single-beam experiment is shown in the third col- 
umn. In this case, the images produce by just one beam (the right 
one) are shown to compare. The result of the double-difference 
described by Eq.|9]is shown in the right column. In the first and 
second columns the polarization pattern is contaminated by the 
intensity term described by 61, and the third column shows that 
the images are more noisy. The fourth column has no artifacts, 
showing a clean polarization pattern, with the lowest noise level. 



6. Additional Corrections 

6.1. Instrumental Polarization 

The polarized light measured by ExPo contains not only the 
polarized light coming from the circumstellar environment, but 
also telescope polarization ( |Miller|[l963| l and interstellar po- 
larization (Hiltner"1949'). Moreover, any small difference be- 
tween the beamsplitter transmission coefficients (Ta,Tb) or the 
FLC transmission coefficients (Tl,Tr) introduces an artificial 
polarization signature. All these contributions are proportional 
to the intensity of the incoming light, so their effects can be 
measured by measuring the degree of polarization at the cen- 
ter of the (largely unpolarized) star The polarization compen- 
sator used by ExPo minimizes these effects, but there is always 
some leftover instrumental polarization. Taking into account that 
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Fig. 9. Uncalibrated Polarized Image. Top rows: AB Aurigae. Bottom rows: Unpolarized Star. Color bar: (uncalibrated) polarization, 
P (Arbitrary Units). The results plotted in the left and center columns are contaminated by instrumental artifacts. These artifacts are 
minimized when applying the double-difference, as is shown in the right column. 



+ ~ Ta ■ 5a * / * 5p, and B]^ + ^ Tb ■ S b * I * Sp, the 
polarization degree at the star's center for the A frame can be 
described as: 



At, + Ab 



cenler 

and for the B frame: 



fii, - Bt! 



5p 



7 



7 



iBias, 



B 



Bias- 



(10) 



(11) 



cemer 

P' refers to any polarization feature that does not come from 



the circumstellar environment. As explained in [Rodenhuis et al. 
( 201 l| l, the polarization compensator is working for all observa- 
tions (including unpolarized diskless stars), and it is impossible 
to set a true zero or reference level for our polarimetric obser- 
vations. Therefore, all calculations are based on the assumption 
that the starlight is unpolarized. However this is not true in case 
there is strong inner disk polarization (Min et al. 201 1 ). In these 
cases, a degree of polarization of the order of a few per cent is 
expected to be measured at the star's position. On the other hand, 
the expected instrumental polarization is a few percent as well, 
so in practice it is impossible to disentangle inner disk polar- 
ization from instrumental polarization. The Asias and Beias are 
then subtracted from the degree of polarization for the A and B 
frames: 



Pa ^ Al-Ar 
/a Al + Ar 



-A 



Bias? 



R 



Pb ^ B^-B 
Ib Bl + Bp_ 



B 



Bias, 



(12) 



(13) 



where dPp, and dPz are the instrumental polarization-corrected 
degree of polarization for the A and B frames, respectively. The 
(uncalibrated) polarization image Pa is then calculated as the 
product of the degree of polarization times the intensity: 



^-(Al+Ar), 



(14) 



Before IP Correcil' 



After IP Correction 









1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 






w 
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Fig. 11. Effect of the instrumental polarization (IP) on the final 
images. Top: AB Aurigae. Bottom: Unpolarized star. Color bar: 
Polarization (P), arbitrary units. 



= 7^ ■ (Bl + Br). 

^B 



(15) 



Figure 1 1 shows the effect of this correction. An artificial bias 
appears in the polarization images when the IP correction is not 
applied. In case the star is unpolarized (bottom row), a bias level 
appears as a consequence of the imperfect compensation. In the 
case of AB Aur (top row) the polarization pattern is distorted 
when no correction is applied. 
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6.2. Field Rotation 



All alt-azimuth telescopes are affected by field rotation. This 
term describes the amount of sky rotation, and it depends on 
the altitude (/i) and azimuth (a) of the observed target, as well as 
the latitude (0) of the observer It will reach its maximum value 
when pointing towards the zenith, and it will be almost negligi- 
ble when pointing to very low altitudes. If the observations are 
made at the Nasmyth focus, an extra rotation must be taken into 
account, due to the relative rotation between the telescope tube 
and the Nasmyth platform (see Marois et al.||2006 Joos||2007 
chapter4): 



G = 



cos ll 



H — — sina ■ cos (p ■ ojq, 
F ^G + H, 



(16) 

(17) 
(18) 



were cjq is the sidereal rate. G represents the rotation rate due 
to the sky movement, and H accounts for the rotation rate due to 
the telescope movement. Equation 18 is the total amount of field 
rotation rate at the Nasmyth platform. The + is related to the two 
platforms at the telescope. In the case of ExPo, the instrument 
has been placed at the "B" platform of the William Herschel 
Telescope, so the "H-" applies. 

While most of the instruments incorporate a de-rotator to re- 
move this effect, the design of ExPo does not include one. It 
is necessary to calculate F and then numerically rotate each ob- 
served target. Due to the short exposure time of the ExPo images, 
the amount of de-rotation between one image and the next image 
is too small to correct (< 0.5°), so rotating each frame individ- 
ually is not necessary. A rotation per image blocks is performed 
instead. F is assumed to be constant within a certain interval. 
If the observed target is not close (< 15°) to the Zenith, this 
interval can be set to 30 seconds. Otherwise, shorter time inter- 
vals must be taken to minimize errors. Different rotation routines 
were tested. To test this, different template images were first ro- 
tated and then de-rotated by using different codes. The standard 
deviation cr of the difference between the original template and 
the de-rotated version is a measurement of the accuracy of the 
rotation procedure. This analysis shows that the best results are 
obtained when the image is first resampled to one fifth of a pixel. 
The rotation is then performed with the IDL routine rot, using 
cubic interpolation, with the interpolation parameter set to -0.7. 

6.3. Sky Bacl<ground Polarization 



The polarization of the sky ( [Bailey et al.||2008l is variable, de- 
pending on the position of the target as well as the position of 
the Moon (since scattered moonlight is polarized). Unlike the in- 
strumental polarization discussed in §6.1, this contribution does 
not depend on the brightness of the observed target, but it be- 
haves like a polarization bias over the whole image. Therefore, 
the background polarization cannot be calculated by measuring 
the degree of polarization at the star's center. The background 
polarization is subtracted at the very last stage of the data anal- 
ysis, by measuring the polarization in the sky regions of the im- 
age. The median value of at least four different sky regions is 
calculated and subtracted from the reduced, polarization images. 

6.4. Ghosts 

Three different sources of optical ghosts are identified in ExPo: 
the filters, the beamsplitter and the FLC. While the latter pro- 



-5 
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Fig. 12. Intensity images of the same star observed without (left) 
and with (right) Hacont filter. The images are plotted at differ- 
ent contrast in order to enhance the ghost image. The diffraction 
caused by the telescope spiders is clearly seen in the left image. 



duces ghosts that move according to the FLC angle, the ghosts 
produced by filters and the beamsplitter remain in a fixed posi- 
tion (an example of these ghosts can be seen in Fig. 12 1. Filter 



ghosts are minimized by tilting the filters. All this ghosts appears 
are contrast ratios of 10"^. 

To analyze the effects caused by the ghosts on the polar- 
ized images, the mean and standard deviation (cr) was calculated 
around the position of the ghosts and around sky (empty) regions 
in the polarization images. This test was performed for several 
(ghost-affected) images, with the result that it is impossible to 
discriminate between ghosts and sky regions when looking at the 
polarized images, i.e., ghosts do not contribute to the polariza- 
tion image any more than the sky background. However, special 
care has to be taken when looking at the degree of polarization 
of the images. In this case, since the ghost appears in intensity 
(at contrast ratios of ~ 10""*), it can contribute to the degree of 
polarization. 



7. Discussion 

Flat-field errors (i.e. pixel to pixel sensitivity variations) are not 
limiting ExPo. Since each image is shifted individually before 
averaging, flat-field effects are averaged out. However, flat-field 
effects might play an important role when an AO system is added 
to the instrument. In this case, the PSF is expected to remain in a 
fixed position on the CCD, with no need to re-center the images. 

Brightest-speckle alignment is not as good as template- 
alignment for our images. This result, in contradiction to pre- 
vious work ( Law|2006 Christou|1991 1 can be explained by the 
characteristics of this instrument. The small field of view used 
by ExPo, combined with the four different PSF's generated in 
each FLC cycle, makes a big difference with respect to the stan- 
dard lucky-imaging observations. The speckle pattern produced 
by a single star is distributed over several pixels, and the beam 
differences make this pattern to be slightly different in the left 
and right images, even though they correspond to simultaneous 
observations. As a result of this, the brightest speckle on the left 
image is not necessarily the same as the brightest speckle on 
the right image. Therefore, the template center method produces 
better results than the speckle center method. However, in case 
of good seeing (below 0.7"), both methods converge, leading to 
very similar results. This can be explained by taking into account 
that, in those cases, the PSF peak is concentrated in very few 
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pixels, so the region covered by the speckle pattern is very small. 
Under these conditions, the center defined by the template-center 
method and the brightest speckle are about the same. 

While theoretically both the ratio and the difference method 
should produce similar results, this analysis shows that the first 
one fails when applied to a dual-beam imaging polarimeter like 
ExPo. The photon noise variance is signal dependent, i.e. it is 
proportional to the amount of photons in each image. On the 
other hand, each image produced by this instrument is affected 
by a different combination of beamsplitter and PLC transmission 
coefficients, which means that the total amount of light in each of 
the four images produced in one PLC cycle is slightly different. 
This also applies to the photon noise, resulting in different levels 
of noise in each image. The double-difference minimizes this 
effect, while the dual-ratio can increase it. This is a particular 
result for this experiment, but it can be generalized to any dual- 
beam polarimeter 

Numerical simulations show that the bias correction de- 
scribed in §6.1 must be applied even when there is no instru- 
mental polarization ( |Min et al. 2011). The strong polarization 
produced by the inner rim of a protoplanetary disk dominates 
the polarized intensity, making it much harder to detect the re- 
maining circumstellar polarization. 

As mentioned in §2, future versions of ExPo may make use 
of 3 PLCs in a Pancharatnam configuration to minimize the 
chromatic effects introduced by the current, single-PLC config- 
uration. This change will not affect the current data reduction 
process, but it will improve the instrument polarimetric perfor- 
mance. 

The inclusion of an AO system in ExPo will lead to a much 
sharper PSP The amount of speckles in these new PSPs will 
be considerably smaller than in our current PSPs. Real speckle- 
sized polarized features (i.e. exoplanets) might be removed as 
a result of the current CR rejection algorithm. One possible so- 
lution to this problem might be a new technique based on the 
comparison of the two differences (AA and AB) described by 
Eq. [5] and Eq. |6] A polarized speckle-sized feature such as an 
exoplanet is expected to appear in both images at the same coor- 
dinates, while in case of a CR, it will appear in only one of these 
differences. 

Chromatic effects introduced by the PLC and the BS are out 
of the scope of this paper, but must be taken into account in 
future instruments such as SPHERE or EPICS. 

Pinally, lucky imaging techniques (i.e. frame selection) does 
not produce significant improvements in our observations, since 
the standard exposure time for a target is around one hour 
(around 20 minutes per PLC position). In this regime, the 
amount of "lucky images" is not enough to produce a good sig- 
nal to noise ratio. 



8. Conclusion 

Dual-beam polarimeters are becoming more and more popular 
among the astronomical community, such as the ones described 
by [Packham et al.| ( |2005[ ), [Masiero et"aL] (|2008) and 'Nagaraju' 
et al. ( 20081 . The high contrast ratios achieved by polarimetry 
makes this technique a promising tool to characterize exoplanets 
as well as circumstellar environments. 

ExPo can currently reach contrast ratios of 10^ at a four- 
meter class telescope, by means of polarimetry. As it is shown in 
this analysis, instrumental artifacts can be minimized by prop- 
erly combining the data produced by a dual-beam instrument. 
The template centering method produces better results than the 



speckle-centering method. The main Umitations of the instru- 
ment are the beamsplitter imperfections and the PLC transmis- 
sion differences, as explained in §5. Puture versions of this in- 
strument will include an Adaptive Optics system as well as a 
coronagraph, which will require new approaches to the data re- 
duction. 
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